Disordered phase of a two-dimensional Heisenberg Model with S* = 1 
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We study an anisotropic version of the Ji — J2 model with S — 1. We find a second order transition 
from a Neel Q = (vr, n) phase to a disordered phase with a spin gap. 



I. INTRODUCTION 

The Heisenberg model can be mapped to the non- linear 
sigma model (NLa) with an additional term due to Berry 
phases. In ID this term has a dramatic effect which 
was first discovered by Haldanei. Integer spin systems 
have a gap in their excitation spectrum, while half-integer 
spin systems are gapless but disordered, with an algebraic 
decay of spin-spin correlations. In the late 1980s it was 
found that Berry phases do not seem to play a role in 
2D for pure Heisenberg HamiltoniansSiilii. Hence, all 
Heisenberg models on the 2D square lattice are on the 
ordered side of the NLa model. This is in agreement with 
a theorem by Dyson, Lieb, and Simon (DLS)^, extended 
to 2D systems by Neves and Peresi, which states that 
for bipartite lattices with S > 1, the ground state is 
ordered. Furthermore, Monte Carlo simulations^ have 
convincingly shown that the S = ^ model has Neel order. 

AfSeck^extented to 2D systems a ID theorem due to 
Lieb, Schultz, and MattisiS (LSM) which applies to half- 
integer spin systems. He argued that since this theorem 
does not apply to integer spin systems, there could be a 
difference between the disordered phases of half-integer 
and integer spin systems in 2D as well. Haldano^ ar- 
gued that if by some mechanism it was possible to drive 
the Heisenberg model into the disordered phase, singu- 
lar topological contributions known as hedgehogs may 
be relevant. Read and Sachde-\»ii carried out a system- 
atic study of the 2D Heisenberg model in the large- 
limit. Their results were in agreement with Affleck and 
Haldane's predictions. They predicted that the nature 
of the disordered phases in 2D is related to the 'spin' 
value. For odd integer spins, the ground state breaks 
the symmetry of rotation of the lattice by 180 degres, for 
even integer 'spins' the ground state does not break any 
lattice symmetry, and for half- integer 'spins', the lattice 
symmetry is broken by 90 degres. 

The usual criticism against such predictions is that 
they are obtained in the limit of large spin and they could 
be invalid for the more relevant cases of small S. It is 
well possible that a different mechanism can emerge for 
small S. The Ji — J2 model is the most popular model in 
which a possible disordered phase has been searched. For 
5 = |, it is generally accepted that for 0.38 < < 0.6, 
this model has a disordered phase. Among the many 
disordered phases that were proposed^-, the columnar 
dimer phase which was predicted in Ref.(— ) seems to 
gain broader acceptance lately. But, we have recently 
arguedi^ that this conclusion, which seems to be sup- 



ported by numerical experiments using exact diagonaliza- 
tion (ED) or series expansions^'', may be incorrect. Large 
scale renormalization group studies on an anisotropic ver- 
sion of the Ji — J2 model show that in the region were 
the disordered phase is expected, physical quantities of 
the 2D model are nearly identical to those of an isolated 
chain. This suggests that there is instead a direct tran- 
sition at ^ = 0.5 between the Neel Q — (tTjTt) and 
Q = (tt, 0) phases. At the transition point, the system is 
disordered with algebraic decay of the correlations along 
the chains and exponential decay in the other direction. 
This state is consistent with the LSM theorem. 



While the case S = ^ has generated numerous 
studies^, other values of S have not been studied to the 
author's knowledge. Thus, the role of topological effects 
in the Ji — J2 model for small S remains unknown. In 
this letter, we propose to study the case S* = 1. We 
will apply the two-step density-matrix renormalization 
groupie (TSDMRG) to study the spatially anisotropic 
Heisenberg Hamiltonian in 2D, 



H = J|| ^ Si^;Si+ij + J_L ^ SijSij+i 

+ Jd ^(SijSi+i_;+i + Si+ijSi^;+i) (1) 



where J|| is the in-chain exchange parameter and is set 
to 1; J± and Jd are respectively the transverse and di- 
agonal interchain exchanges. Although the Hamiltonian 
(|T|) is anisotropic, it retains the basic physics of Ji — J2 
model. In the absence of Jd, the ground state is a Neel 
ordered state with Q = (tt, tt). When Jd ^ J±, another 
Neel state with Q — {tt, 0) becomes the ground state. 
A disordered ground state is expected in the vicinity of 
Jd = In this study, we will only be concerned with 
the transition from Q = (tt, tt) Neel phase to the disor- 
dered phase. The lattice size is fixed to 32 x 33; the trans- 
verse coupling is set to J± = 0.2 and Jd is varied from 
Jd = up to the maximally frustrated point Jd = 0.102, 
i.e., the point where the ground state energy is maxi- 
mal (see Ref.(— )). We use periodic boundary conditions 
(PBC) in the direction of the chains and open boundary 
conditions (OBC) in the transverse direction. This short 
paper will be followed by a more extensive work'^^ where 
a finite size analysis is performed. 
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FIG. 1: Longitudinal spin-spin correlations (full line) and 
their extrapolation (dotted line) for Jd ~ (top), 0.01, 0.02, 
0.03, 0.04, 0.05, 0.06, 0.07, 0.075, 0.08, 0.085, 0.09, 0.101 
(bottom) as function of distance. 



II. METHOD 



We used the TSDMRGi^ to study the Hamiltonian 
(Q. The TSDMRG is an extension two 2D anisotropic 
lattices of the DMRG method of Whitei^. In the first 
step of the method, ED or the usual DMRG method are 
applied to generate a low energy Hamiltonian of an iso- 
lated chain of lenght L keeping mi states. Thus the su- 
perblock size is 9 x mi^ for an = 1 system. Then 
TO2 low-lying states of these superblock states, the cor- 
responding energies, and all the local spin operators are 
kept. These describe the renormalized low energy Hamil- 
tonian of a single chain. They form the starting point of 
the second step in which J± and are switched on. The 
coupled chains are studied again by the DMRG method. 
Like the original DMRG method, the TSDMRG is varia- 
tional. Its convergence depends on mi and m2 , the error 
is given by max{pi, P2), where pi and p2 are the trunca- 
tion errors in the first and second steps respectively. 

Since the TSDMRG starts with an isolated chain, a 
possible criticism of the method is that it could not ef- 
fectively couple the chains. This means that it would 
eventually miss an ordered magnetic phase. The source 
of this criticism is the observation that the DMRG was in- 
troduced to cure the incorrect treatment of the interblock 
coupling in the old RG. But this criticism misses the fact 
that, in the old RG treatment, dividing the lattice into 
blocks and treating the interblock as a perturbation was 
doomed to failure because both the intra and inter-block 
couplings are equal. In the coupled chain problem, how- 
ever, when the interchain coupling is small, it is impera- 
tive as WilsoniS put it long ago to separate the different 
energy scales. 
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FIG. 2: Transverse spin-spin correlations (full line) and their 
extrapolation (dotted line) for Jd = (top), 0.01, 0.02, 0.03, 
0.04, 0.05, 0.06, 0.07, 0.075, 0.08, 0.085, 0.09, 0.101 (bottom) 
as function of distance. 



III. RESULTS 

The low energy Hamiltonian for an isolated chain is 
relatively easy to obtain, we keep mi = 162 states and 
L = 32. For this size the finite size gap is A = 0.4117 
which is very close to its value in the thermodynamic 
limit Ai/ = 0.4105. This is because we used PBC and 
the correlation lenght is about six lattice spacings. The 
truncation error during this first step was pi = 5 x 10^^. 
We then kept m2 = 64 lowest states of the chain to start 
the second step. During the second step, the ground state 
with one of the first excited triplet states with 6*2 = 1 
were targeted. The truncation error during this second 
step varies from /52 = 1- x 10~^ in the magnetic phase to 
P2 = 1. X 10^^ in the disordered phase. This behavior of 
P2 is consistent with previous tests in 5 = i systems in 
Ref.(20) where we find that the accuracy of the TSDMRG 
increases in the highly frustrated regime. 

We have shown in Ref.(i&) for S* = i that: (i) the TS- 
DMRG shows a good agreement with QMC for lattices of 
up to 32 X 33, even if a modest number of states are kept; 
(ii) spin-spin correlations extrapolate to a finite value in 
the thermodynamic limit. However, because of the strong 
quantum fluctuations present for 5* = the extrapo- 
lated quantities were small and thus could be doubted. 
Furthermore, our prediction of a gapless disordered state 
between the two magnetic phases has been regarded with 
a certain skepticism^^-^^-^"^ because it would be expected 
that such a state would be unstable against some relevant 
perturbation at low energies not reached in our simula- 
tion. This is not the case for = 1, where quantum 
fluctuations are weaker. We thus expect larger extrapo- 
lated values Cx=oa and Cy=oo- The results in Fig.J^I for 
the correlation function along the chains, 



Cx — o (^i/2,L/2+lSL/2+a;,L/2-|-l): (2) 




FIG. 3: Order parameter along the chains as function of Jd- 



and in Fig.® for the correlation function in the trans- 
verse direction, 

= 2 (SL/2,L/2+lSL/2,L/2+y>, (3) 

unambiguously show that in the weakly frustrated 
regime, the system is ordered. Despite the strong 
anisotropy, Cx=oo and Cy=oo are not very different. The 
anisotropy is larger for small x and y. But, due to the 
difference in the boundary conditions, Cx seems to reach 
a plateau while Cy bends downward. This behavior of 
Cy is indeed related to the fact that the spins at the edge 
do not feel the bulk mean-field created by other chains. 

As Jd increases, Cx=oo and Cy=oo decreases and vanish 
at Jd^ « 0.085 and Jd, ^ 0.075 respectively. The differ- 
ence in the value of Jd, in the two directions is probably 
due to the difference in the boundary conditions. 

In Fig.Q and Fig.Q we plot the corresponding order 
parameters nix = \/Cx=oo and my = yJCy=aa- The two 
curves display the typical form of a second order phase 
transition. However, we have not extracted any exponent 
because even though we believe our results will remain 
true in the thermodynamic limit, finite size effects may 
nevertheless be important close to the transition. A sys- 
tematic analysis of this region is left for a future study. 

The phase transition is also seen in the spin gap shown 
in Fig.®. In the ordered state, the system is gapless. 
The finite size spin gap is nearly independent of Jd and 
is in the order of the truncation error p2. At the tran- 
sition which occurs at Jd, ~ 0.075, A sharply increases 
and becomes close to the Haldane gap at the maximally 
frustrated point where we find A ~ 0.3854 for Jd ~ 0.102. 

IV. CONCLUSION 

In this letter, we have studied an anisotropic version 
of the Ji - J2 model with S' = 1 using the TSDMRG. 
We find that for a critical value of the frustration, the 
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FIG. 4: Order parameter in the transverse direction as func- 
tion of Jd- 




FIG. 5: Gap as function of Jd- 



Neel ordered phase is destroyed and the system enters 
into a disordered phase with a spin gap. The value of 
the gap at the maximally frustrated point is close to 
that of the Haldane gap of an isolated chain. This dis- 
ordered phase is consistent with the large N prediction 
of Ref.(— ). This study shows the striking difference be- 
tween integer and half-interger spin systems. For S 
the TSDMRG predicted a direct transition between the 
two Neel phases with a disordered gapless state at the 
critical point. Thus, as we have recently foundii, de- 
spite the fact that the mechanism of the destruction of 
the Neel phase is independent of the value of the spins, 
at the transition point, topological effects become im- 
portant leading to the distinction between integer and 
half-integer spins. 
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